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We compute the entanglement entropy of soft gluons in the wave function of a fast moving hadron 
and discuss its basic properties. We also derive the expression for entropy production in a high energy 
hadronic collision within the Color Glass formalism. We show that long range rapidity correlations 
give negative contribution to the production entropy. We calculate the (naturally defined) temper¬ 
ature of the produced system of particles, and show that it is proportional to the average transverse 
momentum of the produced particles. 


I. INTRODUCTION 

High energy hadronic scattering is of great current interest. The structure of the hadronic wave function has 
been the subject of study in the recent years with the aim of better understanding possible manifestations of gluon 
saturation^]. The evolution of this wave function with en erg y, which drives the evolution of high energy scattering 
amplitude is encoded in the so called JIMWLK equation^-lj. 

The stand out property of this Color Glass Condensate (CGC) wave function is that the soft gluons state “grows” 
out of that of the valence modes. More specifically the wave function of the soft modes is determined by the color 
charge density of the valence part of the wave function. This suggests a certain “order” in the system, since the 
soft gluons are entangled with the valence ones. The measure of (dis)order in a quantum system is entropy. In 
particular entanglement entropy measures the degree of entanglement between different subsets of degrees of freedom 
in a quantum state. Entanglement entropy measures both how far the system is from a pure state, but also how close 
it is to a thermal one. It is therefore interesting to ask what is the entanglement entropy of a CGC wave function. 

Although entanglement entropy is a characteristic of a hadronic wave function, it also indirectly carries some 
information about the structure of the final state in collision of this hadron with a hadronic target. For example 
recently we have shown that Bose-Einstein correlations present in an incoming CGC state [i] manifest themselves 
as ridge type correlations between particles produced in the final state of an energetic collision. The CGC based 
calculation of such correlations jb'j provides a possible explanation of the ridge correlations observed by the LHC 
experiments in high multiplicity p-p and p-Pb events|7HlO|. 

In this paper we calculate the entanglement entropy of the soft modes in the CGC wave function. By soft we mean 
the gluon field modes with the longitudinal momentum in the “last” rapidity bin i] < AY ~ l/a s which are still 
energetic enough to participate in scattering at relevant energy. A similar question has been addressed in EH , where 
a definition of dynamical entropy has been proposed. In the present paper we do not use an ad hoc definition, but 
rather directly calculate the standard entanglement entropy of a quantum state. 

We also consider the entropy of the produced system of soft particles. An earlier attempt in this direction is 
presented infill], while a different approach to the entropy production problem can be found in [3, [l4[ (see also the 
review p~5l ] and references therein). The knowledge of the outgoing wave function in the CGC formalism allows us 
to directly calculate the entropy of the final state. We show that the correlations between produced particles give a 
well defined contribution to this entropy. This contribution is negative in accordance with one’s naive expectation 
that stronger correlations mean a more ordered state. We define in a natural way the temperature of the produced 
system of particles. We show that in the weak field limit the temperature is given by T = tt (k±)/2, where ( k ±) is the 
average transverse momentum of the produced particles. 


II. THE DENSITY MATRIX AND RENYI ENTROPY 


At high energy a hadronic wave function has a large soft gluon component. These are the softest gluons in the wave 
function which are energetic enough to scatter on a hadronic target. They occupy the rapidity interval 0 < 77 < AY 
with AY ~ l/a s . In what follows we will calculate the entanglement entropy of this component of the wave function. 
In the CGC approach the hadronic wave function has the form M 

^[a, A] = ip[A]x[a,p\ ( 2 . 1 ) 

where a is the soft gluon field modes, A are the valence modes (with rapidities 77 > AY) and p a (x ) is the color charge 
density due to the valence gluons. As long as the color charge density is not too large a good approximation to the 
soft part of the wave function is a simple coherent state 


x[a,P\ =ex P<( */. 6 “( fc ) [ a « ( fc ) + a a(- fc )]| |0>, ( 2 . 2 ) 

with the Weizsacker-Williams field b’ a (k) = gp a {k) The creation and annihilation operators entering the above 
equation are the gluon operators integrated over rapidity, 

a?(fc) = -4=/ p- aUv,k ). (2.3) 

v AY Jn<AY 27t 

This structure is typical of the Born-Oppenheimer approximation, where the wave function of the fast degrees of 
freedom (soft gluons) is determined by the background of slow degrees of freedom (valence gluons). The valence part 
of the wave function depends on the energy of the process, or in the present context on the total rapidity by which 
the hadron has been boosted from the rest frame, and is subject to JIMWLK evolution^']. At any fixed energy a 
valence observable that depends only on the color charge density 0\p] is calculated as 

(O) = J D[p\W Y \p]0[p] (2.4) 

The rapidity dependence of Wy is determined by JIMWLK evolution. In this paper we will not insist that W[p\ 
solves JIMWLK equation, but instead will consider the (somehwat generalized) McLerran-Venugopalan model fl7j 

W P [p\ = ATe~ fk Pa(k)pai ~ k) (2.5) 

where Af is a normalization factor. 

Note that the soft wave function ea. (12.2D depends on a single longitudinal degree of freedom, i.e. the gluon field 
mode integrated over rapidity in the interval AY. Only this mode is relevant to our discussion in this paper. Rapidity 
dependence becomes significant only when the rapidity interval considered as soft becomes large enough AY > l/a s . 
For such large rapidity intervals additional rapidity dependence of the wave function appears and our discussion would 
have to be amended. We will not consider this in the present paper. 

The reduced density matrix of the soft gluons in the McLerran - Venugopalan model is [s] 

p = Af J D[p] e~ 27 T 3 WfcT p “ (fe)p “ ( — fc) e * /„ |o)(0| 

where we have defined 4> l a (k) = a l a {k) + a^(—k). The integral over the charge density p can be performed with the 
result 




-i <t> 3 


-e 2 


n 

[n^.*j°> <°i^ 




Here we have introduced compact notations: 

<t>i = \a\ a (x) + a°;(x)} ; 


Mij ~ / M 2 (m,u) 


{x-u)i (y- v) 


3 sab 


{x — u ) 2 (y — v) 2 


( 2 . 6 ) 


(2.7) 


Here M bears two polarization, two color, and two coordinate indices, collectively denoted as {ij}. In eci. (12.61) 
summation over discrete and integration over continuous indices is implied. 
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Our goal is to calculate the Von Neumann entropy of the reduced density matrix ea. (12.6l) . As a warm-up exercise 
we first compute the Renyi entropy. To this end we have to evaluate tr[p 2 ]: 


n n 


tr[p 2 ] = £ —<0|e-^^ [ n n <0|^- m & ra , 


m —1 m'=1 


a 4>i 4>j 


10 ) 


( 2 . 8 ) 


Computation of the matrix elements is straightforward as it involves calculation of two vacuum matrix elements, 
which can be performed by summing over all possible Wick contractions. The explicit form of the light front vacuum 
wave function is 


(<£| 0 ) =Ne~^ t 

As this is a Gaussian, it is convenient to introduce 


Gij = <O|e-* M «^&0j|O)/Vdet[7r/(7r + 2 M)\ = [n + 2M] 


-1 

ij 


(2.9) 


( 2 . 10 ) 


where we use the shorthand notation tt = nSijS ab 6 2 (x — y). The matrix elements in ea. (12.81) do not vanish only when 
the total number of </> insertions in each one is even, that is 


n + v! = 2 to 


( 2 . 11 ) 


The result for tr[p 2 ] is a sum of all possible contractions of MG weighted with combinatoric factors. Given that both 
M and G are diagonal in color, the expansion can be organized in terms of color loops: 


L x = tr[MG ]; L 2 = tr[MGMG }; L 4 = tr[MGMGMGMG \.... 

Each trace over color gives a factor N 2 . Performing the averaging over the vacuum state we obtain 


tr[p 2 ] = det 
Using the identity 


7r 

2m .. 

W 1 [ 

_7T + 2M 

n \(2m — n)\ 

m 7i=0 7 


L™(2m - 1)!! + L™- 2 L 4 6(2m - 5)!! + • • • ] 


( 2 . 12 ) 


(2.13) 


2m 

£ 


2m 


1 (2m -!)!! = £ ‘ 


(2m — 1)! 


n\(2m — n)\ 


71=0 V 71 = 0 

we find that the summation over m exponentiates: 


ti!(2m-n)! 2™ 1 (m— 1)! m! 


tr[p 2 ] = exp 


-tr[ ln(l + —)] +2L 2 + 

7r 


~ exp [—2Li + L 2 


After a more careful examination we find that the other terms exponentiate as well 
tr[p 2 } = exp —tr-[ln(l + —+ y j —-— L 2n = exp 


—tr[ln(l + —)] - i*r[ln(l - (2MG) 2 )] 
7r 2 


Finally using the definition ea. (l2.10D we arrive at the closed-form expression 

tr[p 2 ] = exp 

from which the Renyi entropy is found as 


1 

-2^ 


ln(l + ———) 
7r 


<t 2 = — lntr[p 2 ] = £r 


ln(l + —) 
7r 


(2.14) 


(2.15) 


(2.16) 


(2.17) 


(2.18) 
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III. VON NEUMANN ENTROPY 


We are now ready to compute the Von Neumann entropy of the reduced density matrix defined as 

cr E = —tr[p In p] (3.1) 

The following identity is very useful for this purpose 

lnp = lim - (p £ — 1) (3.2) 

e-s -0 e 

We first compute tr[p N ], for arbitrary N, and then take the limit N —> 1 + e. 

Consider the generalization of eq. (12.81) to the calculation of p N . This expression contains a product of N vacuum 
matrix elements of operators that depend on the field <f >. Each one of these matrix elements is calculated independently, 
and thus the fields entering the calculation of different matrix elements can be considered as independent. We thus 
define the multiplet of replica fields <f>f, a = 1,..., TV. Thus we can write 

tr\p N } = + |0) (3.3) 

where now |0) is the light front vacuum of all the replica fields </>“. Notice the nearest neighbor ’’interaction” between 
the replica fields. The replica fields in ea. (13.31) satisfy periodic boundary conditions (j> N+1 = (j) 1 . 

We further rewrite m 



This is a partition function of a spin chain on a replica space lattice. Ea. (13.41) obviously has the discrete translational 
symmetry in the replica space. The ’’action” is diagonalized by Fourier transforming in a: 

1 N N-l 

4> n = — ^ e^ a " </>“; <T (3.5) 

a—1 n =0 

Notice the periodicity relation cj) N ~ n = </>*". The nearest-neighbor interaction in Fourier space reads 



n,m 


Using 

E „— i%3-a(n+m) at X 

e N i d n+mt N 

a. 

we have 

I>“ --0“ +1 ) = ( e ~ i¥n - x ) { e ~ l¥n - x ) = 4iV E sin2 (^ n ) 

an n 


(3.7) 


(3.8) 


tr\p N ] = N n / 2 


det)^ 1 ^ N ^ 2 
2n 


I\i»r 


exp 


-T Eff 


7r + 4 M sin 2 (— n) 


n —0 


N 


where N n ^ 2 is the Jacobian of the transformation (& Using the identity 


7r + 4 M sin 2 (— n) = 

the Gaussian integral in (13.91) is easily computable 


, , 2 tt s 

it + 2M ( 1 — cos(— n) 


tr\p N ] = det^]^ 2 det 


'N-l r 

n 

,n —0 


27T 

7r + 2 M ( 1 — cos(^- n) 


1 - 1 / 2 ' 




(3.9) 


(3.10) 


(3.11) 
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Now we apply (1.396) of Ref. [l8|, with 1 > 2 MG > 0: 

TV—1 r 


n 

n=0 


27T 

7r + 2 M ( 1 — cos(— n ) 

' v N 


N -1 


„ ,27r . 7T 

l- c os(-„) + — 


= (2M) n P 

n =0 L 

= 2 (M) n [cosh(iV arcCosh[l + ^]) 


- 1 


We finally arrive at the following expression for tr[p N ]: 

M 


1 N 

tr{p JSI } = exp<{ - -In2 - — tr 


hn 


1 

— — tr 
2 


In ( cosh(IV arcCosh[l + 


- 1 


The entropy is computed by taking N = 1 + e and keeping the terms linear in e: 

M 


tr[p l+e ] ~ exp < — ^ In 2 — 


In ■ 


1 

~2 tr 


In 


2 M 2 M 


AM 7 r 

1 + — arc Cosh[l + —] 


= 1+ 2 tr 


In — - \ l 

M 


^arcCosh[l + ^] 


Thus we arrive at a closed expression for the entanglement entropy 


cr E = - tr 
2 


M 


4 M 


ln — + \l 1 + — arcCoshjl + —] 


1 


M 


= -tr < ln-h \/1 


4 M 


ln 


1 + 


2 M 


l + \ 1 


(3.12) 


(3.13) 


(3.14) 


4 M\ 


(3.15) 


To understand the basic properties of this expression let us consider translationally invariant case. In this case the 
matrix M is diagonal in momentum space, 


M$(p) = g 2 p 2 (p 2 )^S ab (3.16) 

P* 

In the original MV model /i 2 is a constant and does not depend on momentum, but the momentum dependent 
Gaussian width has been used in recent applications (see e.g. 0)- The contribution to the entropy from large 
transverse momentum modes can be calculated by expanding the expression eg. (13. 151) to leading order in M, since 
for large momenta ( g 2 p 2 < p 2 ) the eigenvalues of M are all small. The expression for the entropy in the weak field 
limit is 


er 


E 

M<1 


tr 


M ne 

— ln — 

7 r M 


(3.17) 


Here we have kept the first two terms in the small M expansion, which are not suppressed by powers of M. 
Thus the dominant ultraviolet contribution is 


’uv 


- -~(N 2 -1)S 

7r 


d 2 p p 2 {p 2 ) 


(27r) 2 


ln 


g 2 d 2 (p 2 ) n( 2 9 2 2/ 2\\ 

6{p - p {p )) 

7r 


e7r p‘ 


(3.18) 


where S is the total area of the projectile. For the original MV model with momentum independent p this expression 
is logarithmically divergent in the UV. Introducing a UV cutoff A we find 


a uv — 7ZZ 2 (JVo 2 — 1)^ 


Aug 2 


In 


2 9 


Q 2 


ln 


9 


Q 2 s 


(3.19) 


where we have identified the saturation momentum in the standard way as Q 2 = g 4 p 2 /ir. The UV divergence in 
this expression is of course the artifact of the eikonal approximation which in the CGC context is applied for all 
transverse momentum modes. In fact the eikonal approximation breaks down when the transverse momentum of the 
soft gluons is of the order of their longitudinal momentum. The natural cutoff on the transverse momentum is of 
order A ~ P + e~ y ~ me Y °, where P + is the total energy of the hadron and m is a typical hadronic scale and Yq is 
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the rapidity of the soft gluons. This cutoff does not depend on the total energy of the hadron, but rather on some 
initial energy at which the eikonal approximation is assume to be valid [38|. 

The contribution from the infrared modes can also be calculated. At small momenta p 2 < Q 2 / g 2 one can formally 
expand in M^ 1 . In this limit the entropy is dominated by the first term in ea. (13.1511 . Here we also keep the first 
subleading correction: 


F 1 r, e2M l 

oo - - Min-] 

z 7r 


and the infrared contribution is 


rfn 


- \(N 2 c -l)S j 


<Pp e 2 g 2 (i 2 (p 2 ) 


(2tt)= 


In 


it p* 


d(Q 2 s - g 2 P 2 ) = ^- 2 (n? - i)SQ 


(3.20) 


(3.21) 


Note that we have chosen to separate the integration region into UV and IR at exactly p 2 = Q 2 /g 2 in eas. (13. 1813. 211) . 
Although the parametric dependence of the separation scale is clear, its exact value is somewhat arbitrary. Our 
reason for choosing the above value, is that at this value of p 2 the integrands in eas. (13.1813.211) exactly coincide. Thus 
this is a unique choice for which the approximation of the total momentum space entropy density by the sum of its 
asymptotic expressions is a continuous function of momentum. 

Combining the two expressions we find approximately for the MV model 


‘uv 


k _SQ 2 2 
aiR -4^ {Nc ^ 1} 


In 


2 9 


Q 2 


In 


? 2 A 2 3 


Q 2 


(3.22) 


A word on the dependence of the entropy on the energy (rapidity) of the hadron. The above calculation is performed 
at fixed rapidity. The hadronic wave function evolves to higher rapidity via JIMWLK evolution. Since this evolution 
is nonlinear, the soft gluon wave function is not given by a Gaussian any more, especially in the saturation regime. 
Calculating the entropy for a nongaussian wave function is significantly more difficult, although we will describe some 
steps toward such a calculation in the last section. However, as has been suggested in the past, approximating the 
soft gluon wave function by a gaussian is not a bad approximation for phenomenological purposes as long as the 
parameters of the Gaussian are taken to evolve with rapidity in a way consistent with JIMWLK evolution[20j. 

In particular one can still use the MV type ansatz but with the function py(p 2 ) taken to be a solution of the 
BFKL[IH (or Bunin) evolution equation. Within this approximation one can derive the evolution of the entropy 
with the total energy of the hadron. In the weak coupling limit, where the evolution of p is given by the BFKL 
equation, differentiating the weak field expression ea. (13. 1 7[) we find 

= -(Wc-l )S J ^2 ln[^(p 2 )]As FKL (p,fc)<?i(fc 2 ) (3.23) 

where cj>(k 2 ) = y M ^ ^ is the gluon unintegrated density, and Kbfkl is the kernel of the BFKL equation. 

In contrast to the MV model, the solution of the BFKL/BK equation at large rapidity exhibits an anomalous 
dimension which slowly varies with transverse momentum 

-PW)~P^ V) \ (3-24) 

7r 

with j(p » Q s ) —► 1; 7 (p ~ Q s ) « .67 and 7 (jp < Q s ) = 0. Note that with such an anomalous dimension the entropy 
diverges as a power of the cutoff in the ultraviolet as a\jv oc A 7 . 

Some properties of the entanglement entropy of CGC are worth noting. First, it is proportional to the transverse 
area of the hadron. This is a natural property of any extensive observable. However cr E is not extensive in the 
longitudinal direction. Recall that a E as calculated above is the entropy on the rapidity slice of width AY, but a E 
is not proportional to AV. The reason for this is of course clear, since only a single rapidity independent mode on 
the interval AV is entangled with the valence degrees of freedom of the hadron. For the same reason the entropy 
does not have any longitudinal UV divergence associated with the fact that we have a sharp boundary between soft 
and valence degrees of freedom in rapidity space. This is very different from a situation one encounters in calculating 
entanglement entropy of a finite region of space A in a local field theory when integrating out degrees of freedom in the 
rest of space B. In the latter case the local interaction between the degrees of freedom along the boundary between 
A and B leads to a UV divergence which depends on the ratio of the correlation length to the (vanishing) width of 
the boundary region. In our case the interaction in the longitudinal direction is nonlocal. The eikonal interaction 
between soft and valence degrees of freedom extends over large distances in rapidity, so that there is no significant 
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contribution to the entropy from the (sharp) boundary. In this sense the entanglement entropy of the CGC is akin 
to the “topological entropy” Q frequently discussed in the context of condensed matter systems, which is a constant 
independent of the boundary or volume of the spatial region under consideration. 

Recently, the “momentum space” entanglement entropy has been discussed in connection with Wilsonian renor¬ 
malization group flow [24]. It would be interesting to explore the relation of these ideas with ea. (13.2311 as a E is indeed 
calculated by separating degrees of freedom in momentum space, while the BFKL equation has been interpreted in 
the past as RG flow in rapidity variable )25j. 


IV. ENTROPY PRODUCTION IN COLLISION. 

Our next quantity of interest is the entropy production during collision of two hadrons. We calculate the entropy 
in the rapidity bin of width AY at rapidity Y away from the valence charges. All the glue at rapidities 77 < Y + AY 
relative to the forward moving particles is considered to be part of the target wave function and as in the standard 
saturation approach is represented by the eikonal scattering matrix S(x). 

According to the eikonal paradigm the effect of the collision is to color rotate the valence as well as the soft gluon 
field in the wave function, so that at time t = 0 right after the collision 


i>[A]x[ a,p] Sip[A]x[a,p] = i/j[SA] x [Sa, Sp\ = ^[5A] exp ji J b l a {x)S(x) [a^{x) + a l a (x)] j |0), (4.1) 

where b°;(x) is the Weszacker-Williams field produced by the rotated color charge density p b (z) = S bc (z)p c (z). Here 
S is the second quantized operator of eikonal 5'-matrix. The action of S to the right (and S t to the left) is equivalent 
to the multiplication of a a (x), a' a (x) and p a (x) by the unitary matrix S(x). Since S is a unitary operator which does 
not mix valence and soft degrees of freedom, it is clear that the initial scattering does not modify the entanglement 
entropy. To see this explicitly we note that in the basis p(x) = S{x)p{x) and (j)(x) = S(x)4>(x) the form of the 
scattered wave function is identical to that of the incoming one in the original basis. Since the integration measure 
for p as well as the vacuum wave function for 4>(x) eq. (12.91) are invariant under this local unitary transformation, 
changing variables to p and </> in the calculation of the entropy immediately establishes that the entropy is the same. 

However the entropy produced in collision is not the same as the entanglement entropy of the final state. Here 
one is interested only in characteristics of inelastically produced gluonic state, and not, say in the contribution of 
soft gluons which are part of the wave function of outgoing “bound state” hadrons. Thus the relevant quantity for 
calculating the production entropy is not the wave function evolved to t —> 00 , but only its “inelastically produced” 
part. This is a similar problem to the one we are faced with when calculating any soft gluon observable produced in 
collision, i.e. single or double gluon inclusive cross section. In that context the solution is well known. Recall, that a 
gluonic observable 0(a,a i ) at t —> 00 is given by the expression [26] 

The multiplication by the extra coherent operator factor accounts for the evolution of the wave function from time 
t = 0 to time t —> 00 as well as restricting Hilbert space for the calculation of O to that of inelastically produced 
gluons only. 

Since this formula is valid for an arbitrary soft observable 0(a, a'), clearly this means that the density matrix that 
describes the produced soft gluons has exactly the same form as ea. (l2.6D but with the matrix M substituted by 

M P = -£5 f u v V 2 (u, v) { ^Zu)2 { (yZv) 2 K 5 M - ^X^M - 5t (y))] ab (4.3) 

Thus the entropy produced in a hadronic collision, which is the Van Neumann entropy of the density matrix specified 
by ea. (l4.3D is given by 


a 


p 



M p 

In-h 

7r 



4 M p , 

-In 

7T 


1 + 


2 M p 



(4.4) 


Although the formal expression for <j p is similar to that for a E , it has a different meaning and its properties are 
significantly different. We stress that a p is not an entanglement entropy like <r E , but rather the entropy of the state 
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produced inelastically in the hadronic collision. As a result a p is UV finite. This is clear, since the UV divergence in 
the calculation of a p came from the coordinate region x —> y —> u —> v, and M p vanishes in this limit, as opposed to 
M. 

Another important property is that a p as defined in ea. (14.4D is calculated for a single scattering event characterized 
by a given eikonal profile S{x). Although it certainly makes sense to talk about event by event entropy production, a 
more global characteristic of scattering would be an event averaged quantity. We will next consider such an average 
production entropy and will relate it to inclusive gluon production amplitudes. 

Note, that eq. (14.31) can be written as a product of single inclusive gluon production amplitude Qi defined in f26j 
averaged over the projectile wave function 

M p = J D[p\ Qi Q 3 W P [p\ = {Qi Qj) P (4.5) 

where 

QUx) = f[ [S(u) S{x)r b p b {u) (4.6) 

2tt J u (x ~ u)- 

with S being in the adjoint representation. 

The target average of this expression is directly related to the single inclusive gluon production probability, or 
rather the phase space density n = dN/dkdydb is 

n{k) = ^ = [ ((QUx)Q?(y)))p,Te ik ^ (4.7) 

d z kdyd z b J xy 

where as usual the averaging over the target fields has to be performed with some weight function W[S']. 

In order to relate the average entropy to inclusive gluon production amplitudes we expand it in powers of the 
fluctuation of the matrix M p around its target average 


M p = M P + (M p - M p )\ 


M p = DSW[S]M P [S ] 


(4.8) 


The target average has the form 


rP _ sab 9 2 [ ..2 r„. (y ~ V )i 


M = — 5 ab / /*>,«) 


47T 2 


{x — u) 2 (y — v) 2 


[Pa{x, y) + Pa{u,v) - Pa{x,v) - Pa(u,2/)] 


where the “adjoint Pomeron” is defined as 


P A {x,y) = 1 - N2 _ i (tr[S t (a:)S(y)])T 


(4.9) 


(4.10) 


Although M p is related to the single inclusive gluon production probability, the two are not equal to each other. 
Assuming translational invariance of /i and Pa we have 


M p = 


S ab 9 2 [ e- ik ^K ij {k,q)p 2 {k-q)P A {q ) 

J k,q 


where 


(k- q)i ( k - q)j h kj (k - q)i kj ( k - q) 3 h 
Pij(P,q) = - - ^771 - 77 + — — - 


(fc - q) 2 ( k - q ) 2 k 2 k 2 (k - q ) 2 k 2 {k - q) 2 k 2 

Conveniently assuming parity invariant target, we can decompose {M p )t in momentum space as 

M P (k ) = U(fc 2 )^ + M t (k 2 ) (Sij - 


k 2 


k 2 


(4.11) 


(4.12) 


(4.13) 


As opposed to M, the matrix M p has both longitudinal and transverse components. These correspond to gluons 
produced with polarizations parallel and perpendicular to their transverse momentum. Clearly both polarizations in 
general are produced, since the gluons in question are free gluons rather than a part of the dressing of the wavefunction 
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of some faster moving valence color charges. The inclusive gluon production cross section [23] is related to the sum 
of the two polarizations 


n{k) oc Mi(k) + M t (k) 

Let us consider the case of small M p , which is relevant for the UV regime. We can then use ea. (l3.17l) 


<t p = (tr 


ire 


M p , 

~ n M p 


)t = tr 


M p ne 
7 r M p 


1 

— — tr 
2n 


|(M P M p ) t - (M p ) 2 |(M i 


Using the representation (14.131) . the first term in (14.151) can be written as 

Mi(k 2 ) 


tr 


M p 7re 

N 2 C ~ 1 

f d 2 k r 

7T n M p _ 

7T J 



Mi{k 2 ) In- 


M t (k 2 ) In 


M t (k 2 Y 


(4.14) 


(4.15) 


(4.16) 


This has the standard form er = — ]Th n* Inn* for the system of noninteracting particle species. The index i here 
refers to transverse momentum as well as the longitudinal and transverse gluon polarizations. Note, that it cannot 
be expressed in terms of the total particle spectrum n(k), unless the longitudinal and transverse eigenvalues of M p 
are proportional to each other. Generically we do not see a reason to expect that this is the case. 

The second term in (14.151) is related to correlated inclusive two gluon production p. [28i - i35[ . Pictorially this term can 
be represented in Fig.la. This diagram (apart from the fact that two of the gluon polarization indices do not close) is 
the same as the connected “glasma graph” contribution to the two gluon correlation [3(i| . The “disconnected” diagram 
in Fig.lb is not included in the expression ea. (l4.15l) . Note that although the diagram in Fig. lb has a disconnected 
topology, it does in fact contribute to the correlated production in the large N c limit if the target averages do not 
factorize^J, 0- I n the language of the Reggeon field theory, this contribution is due to correlated production from 
two Pomerons, while the diagrams in Fig la are due to correlated production from the B-reggeon 37]. 



FIG. 1: 

Two contributions to the inclusive two gluon production. Fig.la - “connected” part that contributes to entropy eg. (14. 151) . 
Fig. lb - “disconnected part that does not contribute to the entropy ea. (l4.15l) . 


With this in mind we can write 


Sa p = — tr 
2n 

j2 I,j2„ 


j(M p M p ) t - (M p ) 2 |(M p )- 1 


(4.17) 


oc / d z kd 2 p {a] a (p)a] D (k)a c j {k)ai{p)) {B ) - (a[ a (p)a“(p))(aJ°(fc)aj(fc)) ( M 0 
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Or in a less convoluted way: 


5<t p = — 


9 


32 7T 5 (TV 2 - 1) 


fj, ( u,v)g (u,v)e 


-\ ik(x-x) 


i,v,u,v,y,x,x,k 

(x - u)i (y - v) m (x - u)j (y - v) r 


kikj 


Mi(k 2 ) k 2 M t (k 2 ) 


u k 2 


(x — u) 2 (y — v) 2 (x — u) 2 (y — v) 2 

x tr([(S(u) - S(x))(S\v) - S\y)MS(u) - S{x)){S*(v) - S\y))]) T + ± [ n{k) 

£ 7T ./u 


(4.18) 


Again, we see that if the longitudinal and transverse parts of M are proportional to each other at all momenta, the 
extra contribution to entropy is proportional to the double inclusive gluon correlation function integrated over the 
gluon momenta. If Mi(k) ^ aM t (k), the longitudinal and transverse gluon polarizations are separately normalized to 
their single inclusive production probability. 

Interestingly the correlation function in eq. (14.171) is normalized by dividing by the first power of the single inclusive 
gluon production and not its square. Recall that the correlated gluon production normalized this way remains finite in 
the CGC approach in the lim it Q 2 S — > oo, while the correlation function normalized to the square of single inclusive 
gluon vanishes in this limit (fl l33l|. 

The natural, but nevertheless interesting property of ea. (l4.15l) . is that the contribution of correlations to entropy is 
negative. This is in accordance with the intuition that a correlated state contains more order and thus has a smaller 
entropy. 

As noted above, the produced entropy is UV finite, since in the UV limit all adjoint Pomerons in ea. (14.91) vanish, 
and the divergent contribution to the momentum integral in ea. (l4.15|) disappears. This is even clearer from the 
momentum space representation of M p ea. (l4.11|) . The perturbative behavior of the adjoint Pomeron at large momenta 
is P(q) oc Qp/q 4 , which leads to M p {k —> oo) — > Q 2 Qp/g 2 k 4 . Even accounting for the BFKL anomalous dimension 
leads to the UV convergent integral in ea. (l4.16l) . The integral in ea. (l4.16l) is also IR finite, since the expected IR 
behavior is P(q —> 0) —> 1 /Q\. 

In general the consequence of eq. (14.911 (neglecting the rotational indices for the moment) is 


" P(?)= *-m F (f) (4 - 19) 

where F is a regular function with the limiting behavior F(x — > 0) —> const x 2 and F(x — > oo) —> const. Using 

eci. (l4.19D we can estimate the entropy production parametrically 


a 


p 


a(N 2 


1)5^| In 

ng 2 


bg 2 Qr 

Qp 


(4.20) 


where a and b are constants of order unity. 


V. DISCUSSION. 

In this paper we have calculated the entanglement entropy of the soft gluon modes in the CGC wave function, as 
well as entropy produced in hadronic collision at high energy. In these calculations we have relied on the MV model 
to represent the valence part of the hadronic wave function. 

The entanglement entropy is an interesting theoretical quantity, as it is a global characteristic of the soft gluon 
density matrix. Its basic properties are determined by the fact that only one longitudinal mode of the soft gluon 
field, namely the rapidity independent mode is entangled with the valence part of the wave function. The entropy 
therefore is not extensive in the length of the rapidity interval, and in this sense is akin to the topological entropy 
discussed these days in condensed matter literature[23[. In a certain sense the analogy is closer than it might seem, 
as the origin of the topological entropy is in the long range quantum entanglement. 

We find that the entanglement entropy of CGC is formally UV divergent and is dominated by the modes with large 
transverse momentum > Q 2 /g 2 . The contribution from the IR modes is finite and suppressed by a square of the 
logarithm of the UV cutoff. 

We have also derived an expression for the entropy production in an hadronic collision. We have shown that it 
is related to inclusive gluon production amplitudes, but is not determined exclusively in terms of these amplitudes. 
In particularly it depends on the production probabilities of longitudinally and transversely (with respect to the 


10 




























direction of their transverse momentum) polarized gluons separately. We also find that correlated gluon emission 
gives a negative contribution to the produced entropy, consistent with the view of entropy as measuring disorder in 
the final state. Just like the entanglement entropy, the produced entropy is not extensive in the rapidity interval over 
which it is calculated. The reason is the same: in the CGC approximation all final state particles are produced from 
a single longitudinal boost invariant mode. 

Having calculated the entropy, the natural question is whether one can define the corresponding temperature. The 
natural definition would be the usual thermodynamic one 


T~ i 


da 

dEl 


(5.1) 


where E± is the transverse energy of the system. 

It is not clear whether this is a sensible quantity for the entanglement entropy. If for E± we take the obvious 

E ± oc f d 2 k\k\M(k ) oc 

J 9 

we obtain for the “entanglement temperature” 

T e oc 


(5.2) 


A 


(JV?-1) In^f 


(5.3) 


which reflects the fact that the entanglement entropy is dominated by the UV modes and is formally UV divergent. 
On the other hand for the produced system of particles ea. (l4.11l) gives 


E ± oc (TV 2 - 1 )S^fQ T (5.4) 

9 

This is just the statement that the average transverse momentum of produced particles is proportional to Qt, while 
the number of particles produced is of the order Q 2 P /g 2 . 

Since the entropy ea. (14.20p is not only a function of E±, the definition of temperature is somewhat ambiguous. To 
define the temperature one has to decide which combination of Qp and Qt has to be kept fixed. Choosing which 
combination to keep fixed corresponds to deciding what is the analog of increasing the energy in the same statistical 
mechanical system. One reasonable choice is to keep Qp fixed, since this corresponds to increasing the energy of the 
produced system by boosting the target, while keeping the projectile unchanged. With this choice using ea. (14.20l) we 
find 


T = cQ T = ^{k±) (5.5) 

where c is a constant of order unity and (k±) = Ej_/N tota i = ^ j ( ppM(l)' > ' s the average transverse momentum per 
produced particle. Note, that even though the relation between (k±) and Qt depends on the details of the function 
F in ea. (14.19[) . the relation between (k±) and temperature is determined exactly. Eg. (15.51) is an intuitively simple 
result showing that the effective temperature is proportional to average momentum per particle. We are nevertheless 
not aware of a direct derivation of this result from the density matrix of the system produced in collision. 

All our calculations were performed in the MV model. The Gaussian weight for the valence charge density averaging 
made it possible to express the entropy directly in terms of a single function M p . The actual weight functional evolves 
with energy and is determined by solving the JIMWLK equation. Although the calculation of entropy in the case of 
non Gaussian W\p\ is more complicated, one can try to approach it in a similar way. 

We would have to calculate 


N 


trlH N ] = tr / J] D\p 0 ]W P \pp] e i ^ k M fc )^( fc )|o)(0|e - *-f* M fe )<£( fe ) 


3=1 


Integrating over pp, and resorting to the same replica trick yields 


/ , , r I \ N/2 r N ( N N 

tr[p N] = ) / n ex P E « + 2 E T ^ a - <^ +1 ] 

^ / J a=1 I z a=1 Z a—1 


(5.6) 


(5.7) 
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Here as before we have introduced N replicas - one for each vacuum matrix elements. In eq. m r is the effective 
action. Interestingly, it is directly related to the projectile scattering matrix, since 


S = / DpWp[p\e 


is L, u p°U) [tf W-0f +1 Wl _ e -|r[0“-0“ +1 ] 


is precisely the scattering matrix element of the projectile on the “target field” — 0“ +1 (u). 

The action T can be expanded in terms of connected correlators of the valence color charge density 


r[0] = r + r % n (\>i 


+ r 4 

1 i. 


(5.8) 


(5.9) 


Tt = 


r?. = 


<5 2 r[</>] , 9- 


5(j)i8(t>j 

8 3 r[4>] 

54>i8(f)j8(t> n 

9 2 r 


\<t >=0 = - 


(x - u)i (y - v)j 


4tt 2 Ju,v ( x - u ) 2 (y - v f 


(p a (u)p b (v)) P 


A=0 = 


9 


( X~u)i (y-v)j ( Z-W) n a b 


(p a {u)p ( v)p c (w)) P 


8tt 3 Ju,v,w (x - U ) 2 (:y - V ) 2 (z - w) 2 

(p“(u)p b (u)p c (u;)/(u/))r n 


(x - u)i ( y - v)j (z - w) n (z' - w')m , at \ b 


l ° nm 16tt 4 J u>VtWtW , (x - u) 2 {y - v) 2 (z' - w) 2 (z' - w') 2 

where the connected correlators are calculated as 

(p 0 (u)p b (x)) P = J D[p\p a (u)p b (v)W P [p\ 

(p a (u)p b (v)p c (w)) P = J D[p]p a (u)p b (v)p c (w)W P [p\ 

(p a (u)p b (v)p C (w)p d (w')) C p mn = [ D[p\ p a (u)p b (v)p C (w)p d (w')W P [p\ \conneoted 


(5.10) 


(5.11) 


If the connected correlators of p are small, these formulae can serve as a starting point for perturbative calculation of 
the entropy. 

An alternative strategy could be to first integrate over <j) in ea. (15.6D . This leads to the following expression: 


N 


N r 


tr\p N ] = / n Dpp exp ^ - J2 


3=1 


3=1 


-In Wp[pp\ 


b 


k.k' 


(bp+i(k) - bp(k)){bp+i(k') - bp{k')) 


(5.12) 


where b^+i = b±. This has the form of discrete lattice model with the same “hopping term” as in Section 2, but with 
a different on site potential V\p\ = — Zn[W [/>]]. 

To take the limit IV —> 1 + e one would need to perform the integration over the replica fields analytically. If the 
projectile is dilute, the potential V\p] is sharply peaked at small p and the integral can be calculated in the “tight 
binding approximation”. For the dense projectile on the other hand, the potential is small since large fluctuations of 
charge density are allowed. In this case one can expand around the free “hopping term”. It would be interesting to 
study the effects of non-Gaussianity on the entropy with the help of these approximations for realistic W|jo]. 

In this paper, in order to compute the entanglement entropy (and the entropy production), we have used the density 
matrix of soft gluons derived after tracing over valence gluons. As a further development of the above ideas, one can 
consider a further reduction of the density matrix. In particular, since gluons with very small transverse momenta 
are unobservable, it could be interesting to trace over gluons with momenta below some soft scale and consider the 
entropy of the remaining modes. Another interesting venue would be to relate two-particle correlations, particularly 
long range rapidity correlations, to the concept of mutual information. The latter is defined through entanglement 
entropy of two different subsets of modes in the hadronic wave-function. 
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